A variety of soliton solutions of time M-fractional: Non-linear models via a unified technique

This work explores diverse novel soliton solutions of two fractional nonlinear models, namely the truncated time M-fractional Chafee-Infante (tM-fCI) and truncated time M-fractional Landau-Ginzburg-Higgs (tM-fLGH) models. The several soliton waves of time M-fractional Chafee-Infante model describe the stability of waves in a dispersive fashion, homogeneous medium and gas diffusion, and the solitary waves of time M-fractional Landau-Ginzburg-Higgs model are used to characterize the drift cyclotron movement for coherent ion-cyclotrons in a geometrically chaotic plasma. A confirmed unified technique exploits soliton solutions of considered fractional models. Under the conditions of the constraint, fruitful solutions are gained and verified with the use of the symbolic software Maple 18. Keeping special values of the constraint, this inquisition achieved kink shape, the collision of kink type and lump wave, the collision of lump and bell type, periodic lump wave, bell shape, some periodic soliton waves for time M-fractional Chafee-Infante and periodic lump, and some diverse periodic and solitary waves for time M-fractional Landau-Ginzburg-Higgs model successfully. The required solutions in this work have many constructive descriptions, and corporal behaviors have been incorporated through some abundant 3D figures with density plots. We compare the m-fractional derivative with the beta fractional derivative and the classical form of these models in two-dimensional plots. Comparisons with others’ results are given likewise.

Fractional derivatives have been widely applied in the analysis of nonlinear evolution equations.It can be implemented to investigate the execution of solutions in various kinds of nonlinear systems.By introducing fractional derivatives, it is possible to express the nonlinear evolution equation in a more general way, which allows for the consideration of nonlinear effects.Fractional derivatives are beneficial in nonlinear models because it can provide more accurate results compared to traditional derivatives, and it hopefully be utilized to analyze solutions with a diverse range of initial conditions.Furthermore, the use of fractional derivatives can reduce the need for tedious parameter changes and simplify the computational process.Recently, many researchers have developed different types of fractional derivatives and applied them to diverse nonlinear systems.M-fractional derivative is used to the Paraxial Wave model [27] and Klein-Gordon model [28], conformable fractional derivative is applied to the date-Jimbo-Kashiwara-Miwa equation [29] and variant Boussinesq equation [30], Caputo fractional derivative is applied to the prey-predator model [31] and Korteweg-de Vries (KdV) models [32], Riemann-Liouville fractional integral is used to the Kudryashov-Sinelshchikov equation [33] and relaxation-oscillation equation [34], fractional beta derivatives is used to Cubic Nonlinear Schro ¨dinger Equation [35] and heat equation [36], etc.
In this paper, we investigate novel truncated M-fractional derivatives via a unified method for the two nonlinear time M-fractional evolution models, such as the Chafee-Infante (tM-fCI) and Landau-Ginburg-Higgs (tM-fLGH) models.There are several methods used for these models in classical differential form [37][38][39][40][41][42][43], but still, the truncated M-fractional derivative on the models are not used.This is the first time of exploration of the Chafee-Infante (tM-fCI) and Landau-Ginburg-Higgs (tM-fLGH) models with such fractional form.The (1+1)-dimensional tM-fCI model equation can be presented as follows: where α is the coefficient of diffusion, which maintains solitary depiction by modifying the balance between the diffusion term and nonlinear component.The C-I equation has some significant phenomena in the homogeneous medium and gas diffusion [38].There are many researchers who have investigated the exact solitary wave solutions through diverse methods such as the exp-function method [37], the direct geometric approach [38], the improved Kudryashov method [39], etc. Besides, a novel nonlinear evolution equation (NLEE) with power law nonlinearity [40] is described as where b 1 ,b 2 ,b 3 ,b 4 are free constants.Due to particular values p = 3,b 4 = 0, the NLEEs [40] reduce to a special nonlinear structure [41], that states Besides this, a typical structure can also be formed from the previous form [41] due to specific parametric constraints b 1 = −1,b 2 = −a 2 ,b 3 = b 2 identified as the LGH model, investigated in [42].
In our view, this model has few gaps to explain accurate real wave nature, and we are going to express more than the above wave, including a truncated M-fractional order derivative on a temporal variable, which can be mentioned as the truncated M-fractional LGH (tM-fLGH) model: where the time and space coordinates are t and x, respectively, and the real constants a and b specify the ion-cyclotron wave and electrostatic potential respectively.The LGH Eq (2) was first developed to characterize the drift cyclotron movement for coherent ion-cyclotrons in a geometrically chaotic plasma.

Preliminary fractional derivative
Characterization: let us consider a mapping H :(0,1)!R, the p order Truncated M-derivative of H exhibit as: Here H d (z) is a one-parameter truncated Mittag-Leffler function that is well-defined as [45,46]:

Unified method
In this section, to explore the unified scheme [27,47], we consider a higher dimensional PDEs: here H = H (x,t).
To solve the Eq (3), the essential steps of the unified scheme are below: Step-01: At first, incorporating the directional wave variable , the PDE (3) becomes an ordinary differential equation.Then we get, Step-02: Consider the solution of Eq (4) is: Where μ S ,ϑ S (p = 0,1,. ..,N) are unknown constants.The trail solution Eq (5) satisfies the Riccati differential equation: To find the balance number N, we use the following formula: Step-03: Injecting different form of Eq (5) into Eq (4), one gets a sequence of } ±h ,h2R.The set of equations is formed if the coefficient of } ±h is thought to be zero.By using Maple 18 to solve the obtained system, the expected solution sets are attained.
Under the conditions on the parameter of Eq (6) the following solution are gotten: For γ<0, ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi À ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi À g p ðZ þ pÞÞ À sinhð2 ffi ffi ffi ffi ffi ffiffi À g p ðZ þ pÞÞ ; : ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi À ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi N.B: For the beta fractional derivative, we used the beta time-fractional derivative is:

Application of UN scheme
In this section, we operate the fractional unified scheme to generate some novel soliton solution for the tM-fCI and the tM-fLGH models.

Time M-fractional Chafee-Infante (tM-fCI) equation
We start the time M-fractional Chafee-Infante equation in the following form: Using truncated M-fractional derivative, and the transformation variable, one ahead to Then the Eq (8) becomes, According to Eq (7), the solution of Eq ( 9) is Here μ 0 ,μ 1 and μ 2 are unfamiliar parameters.Eqs ( 6) and (10) have been placed into Eq (9), and in agreement with step-03, we now possess a system of equations.
To get the required solution set, the above system is solved by using software Maple 18.
Set-01: If we make use of the parameters k ¼ ffi ffi ffi ffi in Eq (10), then we obtain the following solutions.
For γ<0, the hyperbolic solutions are s ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi For γ>0, the trigonometric solutions are s ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi q Gðdþ1Þ l t l and h,ℓ,α,γ,p are arbitrary constants.
Set-02: If we insert the value of the parameters k ¼ ffi ffi ffi ffi (10), then we get the specific solutions.For γ<0, the hyperbolic solutions are ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffiffi q Gðdþ1Þ l t l and h; '; a; g; p are arbitrary constants.
For g > 0; the trigonometric solutions are ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi À ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi q Gðdþ1Þ l t l and h; '; a; g; p are arbitrary constants.

Time M-fractional Landau-Ginburg-Higgs equation (tM-fLGHE)
Let us consider the tM-fLGH model in the following form: Using the transformation variable tðx; tÞ ¼ tðZÞ and Z ¼ kx À o Gðdþ1Þ l t l to Eq (27), we get Let the trial solution of Eq ( 27) is: here m 2 ; m 1 and m 0 are the unfamiliar.In accordance with steps 02 and 03 the Eqs ( 6), ( 28) and (29) provides the succeeding set of equations: To get the required solution set, the above system is solved by using software Maple 18.
Set-01: Inserting the value of the parameters o ¼ 3 2 ffi ffi ffi ffi ffi ffi ffi g p in Eq (10) obtain the solutions as: For g < 0; the hyperbolic solutions are attained, r ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi r À ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi À ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi q Gðdþ1Þ l t l and h,ℓ,a,b,γ,p are arbitrary constants.
For γ>0, the trigonometric solutions are r ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi r À ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi À ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi q Gðdþ1Þ l t l and h,ℓ,a,b,γ,p are arbitrary constants.
Set-02: Substituting the parameters o ¼ 1 2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi we achieve solutions are For γ<0, the hyperbolic solutions are ffi ffi ffi ffi ffi ffiffi À 1 g s ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi 4gk 2 þ2a 2 g q Gðdþ1Þ l t l and h,ℓ,a,b,γ,p are arbitrary constants.
For γ>0, the trigonometric solutions are ffi ffi ffi ffi ffi ffiffi À 1 4g s ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi 4gk 2 þ2a 2 g q Gðdþ1Þ l t l and h,ℓ,a,b,γ,p are arbitrary constants.

Result and discussion
The unified method has been successfully implemented on the truncated M-fractional differential models to obtain some novel soliton solutions.The novel dynamical solutions have been attained by operating this scheme for two powerful models, the tM-fCI and tM-fLGH models.
The parameter of the novel truncated M-fractional derivative has influenced the behavior of attained solutions and also compare the M-fractional derivative with beta fractional derivative and classical form of the differential models.The proposed answers are innovative and have significant implications for illuminating intricate physical nonlinear systems in optical communications.

The time truncated M-fractional Chafee-Infante equation
In this subsector, we graphically demonstrate the impact of fractional parameters on derived soliton solutions of the tM-fCI equation with the three-dimensional graphs and the compare the M-fractional derivative with the beta fractional and classical form of the Chafee-Infante equation in the two-dimensional graphs.The observed solutions express themselves as hyperbolic, rational, and trigonometric function forms with concise descriptions of numerous new types of traveling speeds under varied conditions.The obtained solutions have many significant to describe diverse phenomena in the homogeneous medium and gas diffusion.For γ<0, the unified approach provides hyperbolic function solutions in Eqs ( 11)-( 14).11) through ( 14) provided a variety of structures, such as those seen in Figs 2-4.

The time truncated M-fractional Landau-Ginzburg-Higgs equation
In this subdivision, we graphically illustrate the impact of the fractional parameter [λ = 0.3,0.6,0.9] on the derived soliton solution of the tM-fLGH equation and compare the effect of M-fractional derivative with beta fractional derivative and classical form in two-dimensional graph.Fig 11 to 16 depict the dynamic characteristics of the solutions that were obtained from the tM-fLGH equation.The LGH Eq (2) was first developed to characterize the drift cyclotron movement for coherent ion-cyclotrons in a geometrically chaotic plasma.The obtained solutions have many significant to describe the theoretical framework used in particle physics and

Comparisons and novelty of this manuscript
In this segment compares the attained solutions of Chafee-Infante equations with Habiba et al [39] and Sakthivel et al [37] and solution of Landau-Ginzburg-Higgs equations with Barman et al [44] and Iftikhar et al [46] solutions.The obtained solutions are represented as kink shape and anti-kink shape solutions for the numerical form.Sakthivel et al. [37] discovered the soliton waveform solutions to the classical form of Eq (1) using the exp function technique.They found eight solutions (Please see ref. [37]).Otherwise, we have originated sixteen solutions to Eq (1) by operating the unified method in this article.The extended tanh method and Two variable method for LGH equation Barman et al. [44] explained some important waveforms of the LGH equation using the extended Tanh method and found only hyperbolic solutions.The attained solutions are embodied as bright-type and dark-type soliton, peakon-type, compact, and periodic solutions for the numerical form.Iftikhar et al. [46] discovered the soliton wave solutions to the classical form Eq (2) utilizing the two variable methods and found only two solutions (Please see Ref. [44]).The obtained solutions are represented by two types of singular solutions in the numerical form.Otherwise, through employing the unified strategy in this article, we have discovered sixteen solutions to Eq (2).The Riccati equation is not the same for both methods.
Our Novelty: In this paper, we utilized the unified technique to solve the C-I and LGH models.By using this method, we have found many solutions as trigonometric function, hyperbolic function, and rational function solutions.For the special values of the parameters we get kink shape, the collision of kink type and lump wave, the collision of lump and bell type, periodic lump wave, bell shape, and some periodic soliton waves for tM-fCI model and kink shape, periodic lump wave, and some diverse periodic-and solitary-waves for tM-fLGH.From the above article, it is clearly that some of them are obtained first time for this model.Also at the first time we compared the effect of diverse fractional parameters on the obtained solutions.

Conclusions
In this article, more abundant new exact soliton solutions are successfully developed from two nonlinear truncated M-fractional models, the tM-fCI and tM-fLGH, by applying a unified scheme.The solutions are formed by rational, trigonometric, and hyperbolic functions under the state of kink shape, the collision of kink type and lump wave, the collision of lump and bell type, periodic lump wave, bell shape, and some periodic soliton waves for tM-fCI and kink shape, periodic lump wave, and some diverse periodic-and solitary-waves for tM-fLGH successfully.Exactly, the amplitude and shape of the wave are reformed and changed due to slight changes in the fractional differential order.Moreover, effects of various fractional derivatives are explored in the same 2-D graphics.All the solutions are illustrated with three-dimensional density plots.The output of our research shows that the projected method is an identical, effective, succinct, and strong mathematical tool for integrating complex nonlinear fractional models.In the future, spatio-temporal fractional derivation will be used for these models and also find some novel solitary wave solution by using generalized method.

Fig 1
displays from the numerical form of the derived hyperbolic solutions.In Fig 1, the kink-shaped wave solution is achieved for diverse value of the fractional parameter [λ = 0.1,0.5,0.9] at α = d = 0.5 = −γ, h = 5,ℓ = −4,q=1.The obtained solution is illustrated with three dimensional plots with density.The solutions in Eqs (